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We examine the current-induced magnetoresistance oscillations in high-mobility two-dimensional 
electron systems using the balance-equation scheme for nonlinear magnetotransort. The reported 
analytical expressions for differential magnetoresistivity at high filling factors in the overlapping 
Landau-level regime, which show good agreement with the experimental observation and the nu- 
merical calculation, may be helpful in extracting physical information from experiments. 



In addition to the universally existing Shubnikov-de 
Haas oscillations (SdHO), many different kinds of mag- 
netoresistance oscillations were discovered in the past 
few years in high-mobility two-dimensional (2D) elec- 
tron systems (ES) subject to a weak perpendicular mag- 
netic field and have become a field of great interest. 
These resistance oscillations always occur at low tem- 
peratures and are related to magnetotransport of 2D 
electrons occupying high Landau levels (LLs). Among 
them the microwave-induced magnetoresistance oscilla- 
tions and the related zero-resistance states were the 
central focus of most experimental^^^' 7 -' 8 -'*^ - and 
theoreticalii^i^^^il^i^^ studies. Recently 
the oscillatory behavior in the nonlinear magnetotrans- 
port has attracted much attention: in a 2D system even 
without irradiation, a relatively weak current can induce 
drastic suppression and strong oscillations of the differ- 
ential magnetoresistance, and may result in a state of 
zero-differential resistancej 22 ' 23 i 24 ' 25 i 26 

Several theoretical models have been proposed in 
an attempt to explain this interesting nonlinear 
phenomenon * 27 ' 28 ! 29 ' 30 Numerical examinations based on 
the current-control transport scheme were shown in good 
agreement with the experimental observation for differ- 
ential magnetoresistivity as a function of the ratio of 
the current density to the magnetic field in both the 
magnetic-field sweeping and the current-sweeping config- 
urations covering both separated and overlapping Lan- 
dau level regimes. 27 

It is reported recently that by a systematic analysis 
of current-induced magnetoresistance oscillations, impor- 
tant physical information about electron-electron inter- 
action on the single particle life time can be extracted. 31 
From the point of view of experiment, an analytical ex- 
pression for differential magnetoresistivity, even applies 
only within limited ranges, is highly desirable because it 
can be of great help to extract important physical infor- 
mation from experimental data. So far, a reliable ana- 
lytical expression derived from experimentally confirmed 
theoretical models is still lacking. 

We examine this issue based on the balance-equation 
scheme of nonlinear magnetotransport,— which deals 
with a 2D system consisting of N s electrons in a unit 
area of the x-y plane and subjected to a uniform mag- 
netic field B = (0, 0, B) in the z direction. These elec- 
trons, scattered by randomly distributed impurities and 



by phonons in the lattice, perform an integrative drift 
motion under the influence of a uniform electric field E 
in the x-y plane. For high mobility and high carrier- 
density systems in which effects of electron-impurity and 
electron-phonon scatterings are weak in comparison with 
their internal thermalization, the steady transport state 
at lattice temperature T is described by the electron aver- 
age drift velocity v and an electron temperature T e . They 
satisfy the following force- and energy-balance equations: 

N s eE + N s e(v x B) + f(v) = 0, (1) 
v-f{v)+w(v) = 0. (2) 

Here f(v) — fi(v) + f p (v) is the damping force against 
the electron drift motion due to impurity and phonon 
scatterings respectively, 
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and w(v) is the electron energy-loss rate to the lattice due 
to electron-phonon interactions having an expression ob- 
tained from the right-hand side of Eq. (j4]) by replacing the 
«7|| factor with £l q \, the energy of a wavevector-q phonon 
in branch A. In these equations, n(x) — l/(e x — 1) is the 
Bose function, U(q\\) is the effective impurity potential, 
M (q, A) is the effective electron-phonon coupling matrix 
element, and U2(g\t,u>) is the imaginary part of the elec- 
tron density correlation function at electron temperature 
T c in the presence of the magnetic field. 

The density correlation function n 2 (q\\ , u>) of a 2D elec- 
tron gas in the magnetic field can be written in the Lan- 
dau representation as 

n 2 (g||,^) = ^pE^'( ( B«|/ 2 ) n 2(' 1 ' Ti »' ( 5 ) 

B n,n' 

n 2 (n,n» = ™ f d £ [f(e)-f(s + w)] 

x ImG n (e + uj) ImG„'(e), (6) 

where Ib = yl/jaBj is the magnetic length, C n ^ l+ i(Y) = 
nl[(n + /)!]- 1 r'e~ y [Lj 1 (F)] 2 with L l n (Y) the associate 
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Laguerre polynomial, f(e) — {exp[(e — sf)/T c ] + 
is the Fermi function at electron temperature T e with 
Ef the Fermi level of the system, and ImG n (e) is the 
density-of-states (DOS) function of the broadened LL n. 

The LL broadening depends on impurity, phonon and 
electron-electron scatterings. In a high-mobility GaAs- 
based 2D system the dominant elastic scattering comes 
from impurities or defects in the background or close 
proximity ) 22 i 33 and phonon and electron-electron scatter- 
ings are also not long-ranged because of the screening. 
The correlation lengths of these scattering potentials are 
much smaller than the cyclotron radii R c of electrons in- 
volving in the transport subject to a weak magnetic field, 
but much larger than R c /n for very high (n 3> 1) LLs. 
In this case, the broadening of the LL is expected to be 
a Gaussian formal 

ImG„(e) = -(2 7 r)2r- 1 exp[-2( £ - e n ) 2 /T 2 }. (7) 

In this, e n = (n + ^)u> c is the center of the nth LL, cu c = 
eB/m is the cyclotron frequency with m the effective 
mass of the electron, and T, the half-width of the LL, is 
B 1 / 2 -dependent expressed as^ 

r = (toc/m-.) 1 ' 3 , (8) 

where r s is the single-particle lifetime or quantum scat- 
tering time of the electron in the zero magnetic field. This 
single-particle life time r s , relating to impurity, phonon 
and electron-electron scatterings, is generally tempera- 
ture dependent. 

The DOS function (J7J) for the nth LL is valid and can 
be used in both separated and overlapping LL regimes. 
The total DOS (double spins) for a 2D system of unit 
area in the presence of a magnetic field is given by 
g(e) = — J^ n ImG n (£)/(7r 2 Z 2 3 ), and the Fermi level Ef 
is determined by the electron sheet density N s from the 
equation J d£f{£)g(£) = N s . In the case of high LL fill- 
ing ef is essentially the same as that of a 2D electron 
gas having the same carrier density N s without magnetic 
field, £f = k F /2m (Igf is the Fermi wavevector), and 
v = £f/lo c is the the filling factor. 

When the LL width 2r is larger than the level spacing 
lu c (overlapping LLs), the Dingle factor 

5 = cxp (-it/lo c t s ) (9) 

is smaller than 0.3. Keeping only terms of the lowest 
order in 5 or of the fundamental harmonic oscillation, we 
have the following approximate expression for the DOS: 

g(e) ps — [1 -2dcaa(2ire/w c )]. (10) 

For an isotropic system where the frictional force f(v) 
is in the opposite direction to the drift velocity v, we can 
write f(v) = f(v)v/v. In the Hall configuration with the 
velocity v in the x direction v = (v,0, 0) or the current 
density J x = J = N s ev and J y = 0, Eq. (JTJ yields, at a 
given v, the transverse resistivity R yx = B /N s e, and the 



longitudinal resistivity and differential magnetoresistiv- 
ity as 

Rxx = -f(v)/(N 2 e 2 v), (11) 
r xx = -(df(v)/dv)/(N 2 e 2 ). (12) 

These formulas, expressing the nonlinear resistivity as 
a function of the drift velocity v without invoking electric 
field, are convenient to direct relate to experiments where 
transport measurements are performed by controlling the 
current. This is the basic feature of the balance equation 
approach, 32 in which the drift velocity v of the electron 
system, rather than the electric field, plays as the funda- 
mental physical quantity to affect electron transport. In 
this approach individual (relative) electrons, treated in 
the reference frame moving at velocity v, do not directly 
feel a uniform electric field. The role of the drift velocity 
v is to provide the electron of wavevector q» with an en- 
ergy qu ■ v during its transition from a state to another 
state induced by impurity or phonon scattering. This 
results in an extra frequency gii ■ v in the density cor- 
relation function rL^giucJ + qu ■ v) in Eqs.(3) and (4). 
Since in a magnetic field the density correlation function 
is frequency periodic, II2 (qrn , w) ~ U.2(q\\,u> + u) c ), due to 
periodical LLs, change of the drift velocity v will lead 
to oscillation of related physical quantities. At low tem- 
perature and high LL filling, in view of U2(g\t,u)) func- 
tion sharply peaking around qu ~ 2fcp, the effect of a 
finite velocity v, after the q\\ integration, is equivalent 
to shift the frequency in the II2 function an amount of 
ujj = 2kpv. Thus when ujj varies by a value of uj c , the 
resistivity experiences change of an oscillatory period. It 
is thus convenient using a dimensionless parameter 

_ ujj Imkpv I 8tt m J 
€j = ~ c ~ eB ~\~Nsl?~B ( ' 

to demonstrate the behavior of nonlinear magnetoresis- 
tivity oscillation. 

The numerical examination of impurity-related differ- 
ential resistivity^ shows that r xx = —(dfi(v)/dv)(N 2 e 2 ) 
oscillates with changing ej exhibiting an approximate pe- 
riod Aej ps 1, with maxima at positions somewhat lower 
than integers ej = n (n = 1,2,3, ...) and minima some- 
what lower than half integers ej = n+1/2 (n = 1, 2, 3, ...). 
The predicted results covering both separated and over- 
lapping Landau level regimes, are in good agreement with 
the experimental observation of Ref .[25| in both magnetic- 
field sweeping and current-sweeping cases. 

Simple and accurate analytical expressions of the 
impurity-induced resistivity r xx can be derived for short- 
range potentials in the case of high filling factor v within 
the overlapping LL regime. 

At temperature T <C £f the impurity- induced linear 
(v — » 0) resistivity r xx (0) is given, to the lowest nonzero 
order in S and the fundamental SdHO part, by 

t'I x (0) = R i0 [l + 25 2 - A5D{X) cos(27w)] , (14) 
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where X = 2tt 2 T/uj c , D{X) = X/ sinh(X), and R i0 is the 
low-temperature (T <C Sp) linear resistivity of the 2D 
electron gas in the absence of magnetic field,— which is 
directly related to the transport scattering time r tr or the 
linear mobility (j,q as i?io = m/(iV s e 2 T tr ) = l/(N s efio). 
The expression ([T3)l shows that 7-^(0), though oscillat- 
ing strongly with changing filling factor v, is always pos- 
itive. Both its oscillatory (SdHO) and non-oscillatroy 
parts are r s -dependent through the 5-involved terms in 
the presence of a magnetic field and may change with 
temperature. The linear resistivity i?io of a 2D electron 
system without magnetic field, however, is not affected 
by its single particle lifetime r s and thus temperature 
independent. 

The impurity-related SdHO [the last term in Eq. (TTJ}] 
exhibits minima at integer filling factor v and quickly 
diminishes with rising temperature T ^> ui c /2tt 2 due to 
the temperature-dependent factor D(X). The tempera- 
ture dependence of the non-SdHO part of the impurity- 
induced linear resistivity of the 2D electron system in the 
presence of a magnetic field can appear mainly through 
the temperature change of its single particle lifetime r s 
entering the Dingle-factor coefficient 2S 2 . 

In the case of overlapping LLs and high filling factor is, 
the nonlinear longitudinal differential resistivity r xx can 
be expressed, at temperature T c -C £f, as 



# i0 {l+ 25 2 G(2TT£ j 



(15) 



-46D(X e ) cos(2irv)S{2irej) } , 

in which X e = 2ir 2 T e /uj c , 

S{z) = J (z) - J 2 (z), (16) 
G(z) = J (z) - J 2 (z) - ? - [%h(z) - J 3 (z)] , (17) 

Jk{z) being the Bessel function of order k. With the help 
of power expansions and asymptotic expressions of Bessel 
functions, we have 



S(z)~l-^z 5 



G(z) 
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(18) 
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for z <C 1; and 



S(z) ~ \ — cos 

" 7TZ 



G{z)^\- 

TTZ 
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4/ V 4 



(21) 



for z » 1. 

In the linear limit (ej -> and T c -> T) Eq. ([TS]) re- 
turns to Eq. (H3J). For nonlinear transport Eq. p31) shows 
that the SdHO term, which diminishes when temperature 
T c ^> uj c /2ir 2 , is modulated with an oscillatory factor 
5(27T£j) by the finite current density through the dimen- 
sionless parameter ej. Furthermore, the non-SdHO part 
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FIG. 1: (Color online) G(2-K€j) function (|17|) and its approx- 
imate expression (|19p and (I2ip for small and large argument 

Z = 27T6,-. 




FIG. 2: (Color online) The reduced differential magnetoresis- 
tivity r xx /Rio calculated from Eq. (|15j) (excluding the SdHO 
part) versus e-j = ujj/u c for several fixed values of 7 = 
■k/tsUOj = 0.5, 0.75 and 1. 



of the differential resistivity, which survives the temper- 
ature rising as long as T c <c £f, exhibits strong oscil- 
lation with changing ej, having an oscillation amplitude 
Ar = rl - R i0 as— 



^ = 2<5 2 G(27t £j ). 



(22) 



Since the Dingle-factor coefficient 25 2 varies smoothly 
with changing ej, it is apparent that the oscillatory 
period and the maxima and minima positions of dif- 
ferential resistivity r xx are determined mainly by the 
G(2nej) function. For ej > 0.5, the r 1 oscillates with 

(20) 

tj approximately in a period Ae^ sa 1, having maxima 



around e» 



i-l/8 [n = 
n + 3/8 (n = 1,2,3, 



1, 2, 3, ...) and minima around 
...). Nevertheless, the accurate 
positions of the maxima and minima vary somewhat with 
ej because of the possible change of 25 2 with ej and the 

ey 2 -dependent coefficient at large ej. 

The G(2irej) function pT|) is shown in Fig. 1, together 
with its approximate expressions (| 19[) and (|21[) for small 



and large arguments z — 2nej. The reduced differential 
magnetoresistivity r xx / R\q calculated from Eq. (|15p (ex- 
cluding the SdHO part), is plotted in Fig. 2 as a function 
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of Cj for several fixed values of 7 = (Tr/T s uij) — 0.5, 0.75 
and 1. This differential resistivity derived here for over- 
lapping LLs, turns out to be in good agreement (except 
for the range of tj < 0.5) with more general numeri- 
cal results (covering both separated and overlapping LL 
regimes) of the B-sweeping figure in Ref.[27| and with the 
experimental measurement in Ref.l25l. 

So far numerical results of theoretical model o 27 ' 28 in- 
dicate that the maxima and minima of differential resis- 
tivity are at positions somewhat lower than ej = n and 



Cj = n+1/2 (n = 1,2,3,...) respectively. The present an- 
alytical expression (f2"2"|) with the G(2nej) function ([17]) [or 
its asymptotic expression (|21j) ] derived from the balance- 
equation magnetotransport model, predicting a reliable 
amplitude and phase of the oscillatory differential resis- 
tance for irej 1, is expected to be useful for an accurate 
extraction of temperature-dependent information of the 
single particle lifetime t s or the Dingle factor from the 
current-induced magnetoresistance oscillations in high- 
mobility 2D electron systems at low temperatures^ 
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